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METHOD AND APPARATUS FOR DETERMINATION OF CARDIAC OUTPUT 
FROM THE ARTERIAL PRESSURE PULSE WAVEFORM 

FIELD OF THE INVENTION 
5 This invention relates to the determination of cardiac output (flow 

velocity or flow volume) from the heart through analysis of the pressure or 
diameter waveform in the upper limb (radial, brachial or subclavian artery) or in 
the neck (carotid artery). 
BACKGROUND ART 
10 A clinical and scientific goal for many years has been to measure the 

cardiac output (blood flow ejected by the heart) from arterial pressure. The 
arterial pressure and diameter pulse waveform is created by this ejection, but 
on account of differences in arterial properties with age, differences in the 
pattern of flow ejection of the heart with age and with weakening of left 
15 ventricular muscle, and with change in hear rate under resting conditions, this 
goal has not been attained with precision. 

Additional problems include the effects of wave travel and reflection in 
the arterial tree which cause variable amplification of the pressure pulse 
between central and peripheral arteries. This has been addressed by Kelly and 
20 Fitchett, J Am Coll Cardiol 1992;20:952-63, van Bortel et al., J Hypertens 

2001;19:1037-44, Pauca, Kon and O'Rourke, Br J Anaesth 2004; 92: 651-7, 
and by other methods eg. US patent No. 5,265,01 1 . 

An object of the present invention is to overcome the major problems 
created by aging and cardiac disease, as well as pressure pulse amplification 
25 in the more peripheral arteries, so that the radial or brachial arterial pulse, 

transformed to the aortic pulse, or the carotid or subclavian pulse in the neck, 
can be used to calculate the velocity of blood ejection from the heart into the 
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aorta during each ejection (mean ejection velocity) and over a full beat (mean 
aortic flow velocity). Formulae utilising the known relationship of aortic 
diameter to body height and weight, or direct measurements of aortic 
dimensions by echocardiography can then be used to calculate blood volume 
5 from blood velocity. 

SUMMARY OF THE INVENTION 

According to one aspect of the invention there is provided a method for 
calculating aortic flow velocity from the directly or indirectly measured arterial, 
aortic, or carotid pressure and/or diameter waveform in which the reflected 
1 0 component of the pressure wave is excluded and the peak systolic flow velocity 
V is calculated from the amplitude P1 of the central pressure waveform using 
the Waterhammer formula: 



where C is the aortic pulse wave velocity. 
1 5 Preferably, ascending aortic pulse wave velocity as used in the 

waterhammer formula is measured directly, estimated from the delay from 
wavefoot to first systolic peak or shoulder, or taken from published data, from 
the US National Institute of Aging or another appropriate source for the 
subject's age. 

20 Allowance is made for difference in aortic pulse wave velocity (PWV) 

with arterial pressure so as to apply normalised aortic PWV data to the 
individual. 

Allowance may be made for the reduced aortic velocity in late systole 
with aging, caused by reduced ventricular contractility in late systole, and 
25 attributable to increased left ventricular (LV) load and LV hypertrophy or 
disease. 
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Allowance may be made for the further reduction in aortic velocity in late 
systole, caused by left ventricular weakening and relative change in the heart's 
pumping action from a "flow source" to a "pressure source". 

Average velocity in the aorta is preferably calculated for the period of 
5 ejection and the period of the cardiac cycle, and for other periods of time (eg. 
per second or per minute). 

The preferred method of the invention permits the aortic flow velocity, 
normalised for that individual, to be expressed in terms of volume by multiplying 
by aortic cross-sectional area determined from echocardiography, other 
1 0 methods, of from tables, then to be expressed as volumetric cardiac output per 
minute. 

According to another aspect of the invention there is provided a method 
for determining cardiac output comprising: 

(i) determining the pressure waveform in the ascending aorta, 
1 5 (ii) determining the amplitude (P1 ) of the initial peak of the aortic 

pressure waveform, 

(iii) determining the aortic pulse wave velocity (C), 

(iv) calculating the peak flow velocity (V) using the formula: 

P1 



20 



V = 



1.05*C 



(v) determining the mean systolic flow velocity (Vms) as a 

predetermined percentage of peak flow velocity (V) to allow for 
predetermined factors, 
25 (vi) calculating the mean cycle flow velocity Vmc using the formula: 



Vmc 



= Vms* 



period of systole 
period of cardiac cycle 
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(vii) calculating the cardiac output by multiplying mean cycle velocity 

by the aortic cross-sectional area. 
The pulse wave velocity may be normalised to means arterial pressure. 
5 The mean systolic flow velocity may be reduced to allow for the effective 

heart rate. In a preferred form of the invention the mean systolic flow velocity is 
reduced by 0.9% for each beat per minute above 65 beats per minute. 
BRIEF DESCRIPTION OF THE DRAWINGS 
Fig 1 shows the aortic pressure waveform, 
10 Fig 2 shows the aortic flow waveform over one cardiac cycle, 

Fig 3a to 3f show the waveforms of the various steps of the method of 
calculating cardiac output according to one embodiment of the 
MODES FOR CARRYING OUT THE INVENTION 

In a preferred form of the invention, the pressure waveform in the 
15 ascending aorta is determined by recording the carotid pressure or diameter 
waveform, calibrating this according to the methods of Kelly and Fitchett (J Am 
Coll Cardiol 1992;20:952-63) or van Bortel et al (J Hypertens 2001 ;19:1037-44) 
or similar methods, and taking this as a surrogate of the aortic pressure 
waveform. 

20 Alternatively, a generalised transfer function may be applied to the 

calibrated pressure wave recorded invasively or non-invasively in the brachial 
or radial artery, using the process described in US patent No. 5,265,01 1 or 
other appropriate methods. 

Through use of differentials as described in the above publications or by 

25 other methods, the initial peak or shoulder of the aortic pressure waveform is 
identified and the height of this peak or shoulder (which typically is 90-1 20 
msec after the pressure wavefoot) above the wavefoot itself is calculated. This 
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is taken to represent the pressure wave generated by ventricular ejection 
before the return of significant wave reflection. 

Alternatively, amplitude of this initial pressure peak or shoulder can be 
calculated directly from the radial or brachial pressure waves by exploiting the 
5 known relationship between brachial/ radial and aortic augmentation (Nichols 
and O'Rourke, McDonald's Blood Flow in Arteries, 4 th ed., Arnold, London 
1998; p.368 figure 16.20) and subtracting aortic augmentation from aortic pulse 
pressure. 

The waterhammer formula (V = P/p.C) is used to calculate the velocity 

10 (V) of blood at peak ejection from the pressure wave unaffected by wave 

reflection (P) assuming density of blood (p) = 1 .05. The term C is aortic pulse 
wave velocity and the aortic pulse wave velocity multiplied by blood density is 
aortic characteristic impedance. This can be recorded directly as from delay in 
wavefeet between the carotid and femoral arteries, or indirectly from data 

15 prepared by Lakatta et al from the US National Institute of Aging (Lakatta EG. 
Cardiovascular Aging in Health. Heart Failure in the Elderly. In: Clinics in 
Geriatric Medicine 2000;16:419-43), which link aorta pulse wave velocity with 
age, and show no gender difference. 

Other normal values e.g. Avolio etah, Circulation 1983;68:50-58; or 

20 Nichols and O'Rourke 1998, may be used instead. Examples of these data are 
shown in Fig 3. [For the Lakatta data, the following formula is used:- C = 8.52 * 
age + 222, where C is aortic pulse wave velocity in cm/sec]. Pulse wave 
velocity is normalised to mean arterial pressure of 100 mmHg by formula:- C 
adj = C - 7.1*(100 - mean pressure), derived from the known passive decrease 

25 in aortic pulse wave velocity with fall in mean pressure described by Asmar et 
al, Hypertension 2001;38:921-6. 
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The duration of ejection from wavefoot to cardiac incisura is measured 
and compared to the total cycle length. This allows for differences in the 
duration of left ventricular contraction and relaxation as seen under different 
conditions eg. change in heart rate. The incisura is determined from the 
5 recorded radial, brachial, subclavian or carotid waveform, using differentials or 
other methods. The ejection duration is typically in the region of 250-350 msec, 
and relative duration of systole 30-40% of the cardiac cycle. 

Assuming no flow in diastole, one then has to consider the shape of the 
aortic flow wave, and how this is influenced by aging and weakening of the 

1 0 ventricular muscle. In the first aspect, it is assumed that peak flow is achieved 
at the already identified shoulder and that flow ceases at the identified incisura. 
It is also assumed that the area under this curve in a healthy normal person 
under age 60 equates to a rectangle characterised horizontally by the duration 
of ejection and vertically by 80% of peak forward flow velocity. It is further 

1 5 assumed no flow in diastole. This step is necessary on account of 
intermittency of cardiac contraction and relaxation. 

The effects of aging are allowed for by assuming relatively lower forward 
flow velocity in the latter part of systole after the early flow peak, as shown by 
Nichols et al (Am J Cardiol 1985;55:1 179-84). As a first approximation, the 

20 value of 80% is taken to reduce by an absolute value of 1 0% for each decade 
over age 60 - ie. to 60% at age 80 and 50% at age 90. More precise 
approximations can be made when more normal human aging data becomes 
available. 

The effect of impaired left ventricular contraction on the pulse waveform 
25 (Westerhof and O'Rourke. J Hypertension 1995;13:943-52) is allowed for by 
reducing the value for age calculated above by absolute 10% if LV ejection 
fraction (LVEF) is known to be between 25-40% (ie. from 80% at age 60 to 
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70%) and by absolute 20% if LV ejection fraction is known to be below 25% (ie. 
from 80% at age 60 to 60%). If cardiac failure is present and ejection fraction is 
not known, the first figure of 1 0% (as for LVEF 25-40%) is used. 

Allowance can be made for the effective heart rate on flow waveform 
5 by reducing the predicted value of mean systolic velocity on mean cycle 

velocity by 0.9% for each beat per minute above 65 beats per minute. This 
correction is applied to the value already corrected for age and the presence of 
cardiac failure. 

Calculated flow velocity is thus determined per ejection (i.e. per beat). 
1 0 This is normalised for body size since the waterhammer formula relates 
pressure change to velocity change. In each individual, velocity can be 
converted to volume flow (stroke volume in mis/beat) from measurement of 
aortic cross-sectional diameter and area by ultrasound, or from the nomogram 
from Lakatta et al (Lakatta EG. Cardiovascular Aging in Health. Heart Failure in 
15 the Elderly. In: Clinics in Geriatric Medicine 2000;16:419-43) which relates body 
height and weight to aortic cross-sectional area. 
From the Lakatta data, the formula:- 
D = (0.0654* age) + 12.63 
can be used where D is diameter per square meter body surface area. 

20 

Body surface area (BSA) may be calculated from the equation of Du 
Bois and Du Bois:- 

BSA = Mass 0 426 (Kg) * Height 0 725 (cm) * 71 .84. 

See Dubois and Dubois (Arch.) Intern. Med 1916; 17:863) or as 
25 converted into nomograms in the Geigy Scientific Tables (-see Fig. 5). 

Cardiac output may be calculated as mean aortic flow velocity times 
calculated aortic cross-sectional area. On account of difficulties in accurately 
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measuring or estimating aortic cross-sectional area, a preference is given for 
expression of average flow as velocity over the whole cycle (this is the same as 
average velocity flow over a full minute). 

The method described has other aspects, and it will be appreciated that 
5 variations and additions are possible within the spirit and scope of the 
invention. Previous applications of this type of method have used the 
waterhammer formula on central or peripheral pressure pulse waveforms. The 
specific novelty of this approach is that ft takes account of the change in left 
ventricular ejection pattern which occurs with aging (as a consequence of early 

10 wave reflection and late systolic pressure augmentation), and the progressive 
change of the heart from a "flow source" to a "pressure source" with aging as 
afterload increases (Nichols et al. Am J Cardiol 1985;55:1179-84), and other 
further change which occurs as the ventricular muscle weakens through 
ventricular hypertrophy or from intrinsic cardiac (eg. coronary) disease 

15 (O'Rourke MF. Blood Pressure 1994;3:33-37). The approach also concentrates 
on scaled values of peak and mean flow in terms of linear velocity rather than 
volume flow. 

Other aspects of the invention allow direct measurement of aortic pulse 
wave velocity, or its estimation from the time to return of the reflected wave 

20 (London et al. Hypertension 1992;20:10-19), instead of from the NIA or other 
tables. Other aspects entail use of the calibrated carotid or upper limb 
pressure pulse waveform to determine the impulse generated by ejection 
(without reflection) and of ultrasonic or other techniques to calculate aortic 
cross-sectional area. On account of different systematic effect of age on the 

25 proximal aorta that on the truncal arteries, it may be necessary to introduce a 
scaling term in the equations used to calculate peak velocity from pressure 
under different conditions. It may also be necessary to vary the constants 
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presently proposed to describe the effect of aging and of heart failure on the 
arterial pulse. 

It is desirable to allow for change in the flow waveform contour at 
different heart rates at rest by reducing the calculated average flow velocity 
5 below that initially calculated when heart rate increases above 65 beats/minute. 
Measurement of cardiac output or cardiac index (scaled to body 
dimensions as achieved here) is an important clinical measure, being increased 
in exercise and pregnancy, with thyroid overactivity, and with some forms of 
"hypertension", including "white coat hypertension". Cardiac index is reduced 
10 with blood or fluid loss, with pulmonary embolism and cardiac failure from 
multiple causes. Since the major function of the heart is to pump blood, 
measurement of the amount pumped is a valuable clinical sign. To determine 
this simply with logical physiological processes is an important advance. 

Fig 1 shows a central (aortic or carotid) pressure waveform recorded 
15 directly or synthesised from a peripheral pressure waveform. The point 0 is the 
foot of the wave from which pressure rises smoothly up to a localised peak or 
shoulder at point 1 , some 90-1 20 msec after the foot of the wave. Pressure 
may rise further after this point 1 , but then declines to an inflection or incisura 
typically some 250-350 msec after the wavefoot. This incisura denotes closure 
20 of the aortic valve and the end of ventricular ejection. 

The pressure rise from point 0 to point 1 (P1 - P0) is determined by the 
flow ejection in the time T1 -TO, by aortic pulse wave velocity and blood density 
according to the waterhammer formula. The flow velocity change from zero at 
the wavefoot to its peak and back to zero at the incisura is represented by the 
25 broken line. The peaks of pressure and flow correspond at point 1 . 

The aortic flow velocity waveform shown in Fig 2 in systole normally 
approximates to a rectangle (dotted lines) of base corresponding to ejection 
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duration and height corresponding to 80% of peak flow. This alters with aging 
as left ventricular load increases and the ventricle hypertrophies. Ejection 
duration may increase, but the patterns of late systolic flow changes such that 
there is low velocity in later systole and average systolic flow typically falls to 
5 60% of peak flow. 

In cardiac failure the same phenomenon is seen with lower flow in late 
systole. The period of systole may also be decreased but this is directly 
measurable as reduction in ejection duration, and of ejection duration/ duration 
of cardiac cycle. Presently, the subtle change in the flow pattern can only be 
1 0 inferred from the lesser degree of augmentation than expected in late systole. 

A nomogram of the relationship between aortic pulse wave velocity and 
age is set forth in NIA studies (Lakatta EG. Cardiovascular Aging in Health. 
Heart Failure in the Elderly. In: Clinics in Geriatric Medicine 2000;16:419-43). 

A nomogram of the relationship between aortic diameter and body 
1 5 surface area is set forth in NIA studies (Lakatta EG. Cardiovascular Aging in 
Health. Heart Failure in the Elderly. In: Clinics in Geriatric Medicine 
2000;16:419-43). 

A nomogram for the calculation of body surface area is included in the 
Geigy Scientific Tables. 
20 The steps of one method of the invention are as follows with reference to 

Figs. 3a to 3f:- 

Step 1 : calculation of pressure rise in aorta caused by peak flow 
ejection (P1) (multiple methods possible), - Fig. 3a. 

Step 2: calculation of peak flow velocity (F) corresponding to (P1) 
25 using water hammer formula. Determination of aortic pulse 

wave velocity from normative data, and correcting this to a 
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value appropriate to the individual's mean arterial pressure, - 
Fig. 3b. 

Step 3: calculation of mean flow velocity in systole, by assuming that 
mean systolic flow = 80% of peak flow - Fig 3c. 

Step 4: allowance for effect of age on flow pattern, with reduced flow 
in late systole causing mean systolic flow to be < 80% of 
peak, Fig. 3d. 

Step 5: allowance for effect of heart failure, with flow in late systole 

further reduced at any given age to « 80% of peak, - Fig. 3e. 

Step 6: allowance for heart rate by further reduction in systolic flow to 
«<80% with heart rate over 65/ minute. 

Step 7: calculation of mean cycle flow velocity as (mean systolic 
velocity) x Z , - Fig. 3f. 



15 



where Z = 



period of systole 
period of cycle 



20 



Step 8: 

Example 

Subject: 



25 



calculation of volume flow from mean cycle velocity and aortic 
cross sectional area as cardiac output. 

Male, 67 years, height 168cm, weight 76Kg. 

•Mean blood pressure: 87mm Hg 

•Initial peak amplitude (P1): 33mm Hg 

•Heart rate: 67 bpm 

•Period of systole/period 36% 
of cardiac cycle: 
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Utilizing the steps set forth above, the cardiac output in litres/minute is 

calculated from the initial peak amplitude (P1) as follows:- 

(i) Pulse wave velocity (C) = 8.52* Age + 222 

8.52 * 67 + 222 
= 792.7 cm/sec 



10 



(ii) Pulse wave velocity 

normalised to mean 
blood pressure 



C-7.1 (100-mbp) 

792.7-7.1 (100-87) 
700.4 cm/sec 



15 



(iii) Peak flow velocity (V) 



P1 



1.05 *C 
33(980*1.36) 
1 .05 * 700.4 
59.8 cm/sec 



20 



25 



(iv) Age adjustment factor = 
(normal cardiac function) 
(AAF) 



(v) Mean systolic flow velocity = 



80% minus 1% for 
each year over 60 

80% - 7% 
73% 

73% Peak flow velocity 
0.73 * 59.8 
43.65 cm/sec 



30 



(Heart rate -65) 90% 
(67-65) 0.9 
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1.8% 



(vii) Adjusted mean systolic 
flow velocity (Vms) 



(100%-1.8%)*43.65 
42.86 cm/sec 



10 



(viii) Mean cycle flow 
velocity (Vmc) 



Vms * period of systole 

period of cardiac cycle 

42.86 * 36% 
15.43 cm/sec 



15 



20 



(ix) Aortic cross-sectional 

area (A c/s) 

(x) Cardiac output 



7.8 

Vmc * A c/s 
15.43*7.8*60 
120.3ml/sec*60 
7.22 Litres/minute 



Various modifications may be made in details of the method without 
departing from the scope and ambit of the invention. 



